Density of states in a superconductor carrying a supercurrent 
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We have measured the tunneling density of states (DOS) in a superconductor carrying a su- 
percurrent or exposed to an external magnetic field. The pair correlations are weakened by the 
supercurrent, leading to a modification of the DOS and to a reduction of the gap. As predicted by 
the theory of superconductivity in diffusive metals, we find that this effect is similar to that of an 
external magnetic field. 

PACS numbers: PACS numbers: 74.78.Na, 74.20.Fg, 74.25. Sv 



How is the superconducting order modified by a su- 
percurrent? The superconducting order is based on pair- 
ing electronic states which transform into one another 
by time reversal 0] . The ground-state wavefunction cor- 
responds to a coherent superposition of doubly empty 
and doubly occupied time- reversed states, in an energy 
range around the Fermi level given by the BCS gap en- 
ergy. When an external magnetic field B — curl A is 
applied, time-reversed states are dephased differently, re- 
sulting in a weakening of superconductivity. In presence 
of a supercurrent, the superconducting order no longer 
corresponds to the pairing of time-reversed states, which 
results in a kinetic energy cost, and again in a weakening 
of superconductivity. In the early stages of the theory of 
superconductivity, it was found that, in diffusive super- 
conductors (in which the electron mcan-frec-path is short 
compared to the BCS coherence length) and in homoge- 
neous situations, the modification of the superconducting 
order by a magnetic field, by a current and by paramag- 
netic impurities could be described by a single parame- 
ter, the depairing energy F (2L_ Later on, the reformula- 
tion of the theory by Usadel Q in the diffusive limit 
extended this equivalence to inhomogeneous situations, 
where the modulus of the order parameter may vary in 
space. In the Usadel equations, all physical quantities 
only involve the intrinsic combination V(/3 — {2e/h)A, 
where the gradient Vtp in the phase of the superconduct- 
ing order parameter is associated with the supercurrent, 
revealing the equivalence of a supercurrent and of an ap- 
plied magnetic field. The Usadel equations are now at 
the basis of the understanding of mesoscopic supercon- 
ductivity in diffusive conductors |^, Q . Experimentally, 
measurements of the density of states (DOS) in a thin 
superconductor placed in an in-plane magnetic field were 
well accounted for by the concept of depairing energy 
0. In contrast, the effect of a supercurrent has been 
partly addressed in a single experiment, focused on the 
reduction of the superconducting gap close to the critical 
temperature . A complication of the experiments with 
a supercurrent is that, if the sample width exceeds the 
London penetration length A^, the current distribution 
given by the non-local equations of electrodynamics 9] is 



not homogeneous. In the experiment reported here, the 
superconductor is wire-shaped, with thickness and width 
smaller than Ai,, so that the current flow is homogeneous 
and the magnetic field penetrates completely. Moreover, 
the effect of the magnetic field induced by the supercur- 
rent is then negligible. This simple geometry allows to 
test the fundamental equivalence between the effect of a 
magnetic field and of a supercurrent in a diffusive super- 
conductor, and to compare precisely with the predictions 
of the Usadel equations. 

Our experiment was performed on a current-biased su- 
perconducting wire made of aluminum, placed in a per- 
pendicular magnetic field B (see Fig. 1). The density 
of states in the wire was infered from the differential 
conductance dl/dV {V) of a tunnel junction formed be- 
tween a small section of the wire and a normal probe 
electrode made of copper. Disregarding Coulomb block- 
ade and temperature effects (see below), dI/dV{V) is 
proportional to the DOS n{eV). The sample was fabri- 
cated in an electron-beam evaporator in a single pump- 
down, using three-angle shadow-mask technique through 
a PMMA suspended mask patterned using e-beam lithog- 
raphy 10\. The substrate was thermally oxidized silicon. 
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FIG. 1: Layout of the experiment: a 10 /^m-long, 120 nm- 
wide and 40 nm-thick superconducting (aluminum) wire can 
be current biased at Is or exposed to a magnetic field B. A 
normal probe electrode forms a tunnel junction (dashed area) 
with the wire. To a good approximation (see text), the differ- 
ential conductance of the junction dl /dV{V) is proportional 
to the DOS in the superconductor. 
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The 10 /xm-long aluminum wire, with width w = 120 nm 
and thickness t = AO nm, was superficially oxidized in 
order to form a tunnel barrier with the copper probe 
electrode overlapping it on an area 150 x 60 nm^. The 
sample was mounted in a copper box thermally anchored 
to the mixing chamber of a dilution refrigerator. Mea- 
surements were performed at 25 mK. Electrical connec- 
tions were made through filtered coaxial lines. From 
the low-temperature, high-magnetic-field wire resistance 
in the normal state, R — 77 fl, the conductivity a = 
27 ri~^./im~^ is infered assuming that the electrical cross- 
section of the wire is S = wt. The diffusion coefhcient 
D = 4:9 cm^.s"^ is then deduced using Einstein's relation 
cr = N{0)e^D, where N{0) = 2.15 lO^^ J-i.m'^ is the 
density of states at the Fermi level of aluminum in its nor- 
mal state and e the electronic charge. The superconduct- 
ing gap Aq — 205 ^eV was deduced from the differential 
conductance- voltage characteristic dI/dV{V) measured 
at S = 0, /s = (dashed line in Fig. 2). Using these 
parameters, we obtain the superconducting coherence 
length ^0 = \/fiD/Ao ~ 125 nm and the London length 
Xl = ^/h/~(|M3^TaAo) ~ 175 nm. Since Xl ^ w/2, the 
current density is homogeneous when the wire is current- 
biased, and a magnetic field penetrates uniformly in the 
wire. The measured critical current of the wire at i? = 
was Ic = 106 /iA. 

In Fig. 2, two dI/dV{V) curves are shown, respec- 
tively measured at Is — 70 /iA, zero field, and at zero 
current, B — 23 mT. The reduction of the gap and the 
smearing of the peak near the gap energy are similar 
in the two situations, bringing already evidence of the 
equivalent effect of Is and B. Note that the magnetic 
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FIG. 2: Measured dI/dV{V) for different combinations of 
the bias current and magnetic field: dashed line: Is ~ and 
5 = 0; solid lines: Is = 70 fiA and B — 0, and Jg = and 
B = 23 mT. 



field created by the supercurrent has a negligible effect: 
for Is = 70 iiA in the wire (see Fig. 2), ^ - 0.15 mT 
whereas the resulting DOS is recovered at Ig = with 
B = 23 mT. A complete set of data is presented in Fig. 3, 
with dI/dV{V) measured for Is = 17, 51 and 85 /iA at 
5 = 0, and for B = 11.5 to 69 mT by steps of 11.5 mT, 
at Is = 0. Note that when the wire is current biased, the 
superconducting state is metastable. In practice, for bias 
currents larger 85 /iA, the system switches to the resis- 
tive state during the recording of the dI/dV{V) curve. 
The measured curve is then similar to that obtained in 
the normal state. In order to account quantitatively for 
the data, we use the Usadel theory 0, Q . In this theory, 
correlations between electrons of opposite spins and mo- 
menta are described by a complex function 9{f,E), the 
pairing angle, which depends both on space and energy, 
and a local complex phase (p{ f, E). The local density of 
states is given by n{r,E) = iV(0)Re[cos(6'(f, £;))]. The 
pairing angle and the complex phase obey the Usadel 
equations: 



^V'^e+[iE-^ ( V - ^1)2 cos 6] sin + A cos = 



V[(V(^- yi*)sin2 



6] = 0. 



(1) 



(2) 



A term describing spin-flip scattering, which is found neg- 
ligible in our experiment, has been omitted here. The 
pairing potential A(r) is determined self-consistently by 



A(r) = 7V(0)V;// 



d£'tanh( 



E 



2kBT 



)Im(sin6l) (3) 



where 14// is the pairing interaction strength, ojd the 
Debye pulsation, ks the Boltzmann constant and T the 
temperature of the superconductor. 
The supercurrent density j is given by: 



' E ^ 2e ^ 

dE tanh(— )Im{sin^ e){V if - —A). 

2kBT h 

(4) 

In a situation like ours where the system consists entirely 
of a single superconductor, V(/3 does not depend on en- 
ergy, and j can be written as a product of the density of 
charge in the superconducting state ps{r) — eN{0)Us{'r), 
with Usir) = d-Etanh( y,.^^ )Im(sin^ 6), and of a su- 

perfluid velocity Vs 



-)Im(sin^ 6 
{2e/h)A). 



■■2kBT) 

D{Vip- 

We have first checked numerically that the dependence 
of 6 on the directions transverse to the wire could be ne- 
glected because the width and thickness are smaller than 
the superconducting coherence length ^q, which is the 
characteristic lengthscale for the variations of 9. As a 
consequence, all the quantities can be replaced by their 
values averaged on the transverse directions. In the Lon- 
don gauge, the effect of the magnetic field is described 
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covers the generic equation given in Ref. 0: 
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FIG. 3: Normalized differential conductance dl /dV{V) of the 
probe tunnel junction: Top: at _B = 0, as a function of the 
supercurrent 7s (from right to left: 7s — 17 nA, 51 fiA, and 
85 /iA). Bottom: at 7s = 0, as a function of the magnetic 
field B (from 11.5 mT to 69 mT by steps of 11.5 mT). Solid 
lines are best fits with dl /dV{V) calculated with an electronic 
temperature dependent on V (see text); dashed lines are the 
best fits with dl /dV{V) calculated with a constant electronic 
temperature. Insets: depairing energy r(in units of the gap 
Ao at 73 = and 7s = 0) for different currents and magnetic 
fields, deduced from the fits of dI/dV{V). In the top inset, 
square symbols correspond to the data in the main panel (73 = 
0), whereas triangles and disks were obtained from data taken 
at B = 10.2 mT and B = 27 mT, respectively. Solid lines: 
Fits with theory, leading to depairing current and magnetic 
field 7r = 240 fiA and Br = 105 mT. 



by a vector potential parallel to the wire axis x, with an 
amplitude A^: = By, so that {A^} = and {A"^) = 
Bw/{2\/3) 0|. The constant phase gradient dip/dx is 
given by the supercurrent Is = jS = UsL/ (eR) {dip/dx). 
Since d'^(p/dx^ = 0, Eq. ^ reduces to i9(sin^ 0)/dx = 0. 
No spatial dependence remains in Eq. (Q, and one re- 



£; + ircos6' = iA 



cos 9 



smf 



(5) 



where 



is the depairing energy, which contains the effect of both 
a phase gradient and a magnetic field. Note that since 
r/Ao = ^{^od^/dx)^ + ^{CowB/{h/e))^ the relevant pa- 
rameters are the phase difference between two points of 
the wire distant by said the number of flux quanta 
in an area wS.o- The depairing energy is related to the 
external parameters 7^ and B by the equation: 



r 



/ Ao 7, 



B 
B^ 



where we have introduced the characteristic depairing 
current and magnetic field 7r = \/2Ao/{eR{£,o)), with 
7?(^o) = R^o/L the resistance of the wire on a length 
and Br = V6(/i/e)/(«;^o)- 

Since the transverse dimen- 
sions of the wire are smaller than the London length Ai, 
the depairing energy due to the induced field is negligi- 
ble (smaller by a factor ~ 10"''' ^2]) compared to the one 
due to the supercurrent. The DOS for a given depairing 
energy F is obtained from the self-consistent solution of 
Eqs. lO and Q. For practical purposes, we give the 
approximate expressions for the resulting A(F)/Ao and 
C/,(F)/Ao, valid, at keT < A, for F/Aq < 0.3: 



A(F) 
Ao 

Ao 



1 - 0.75 



tt/2- 1, 



Ao 
Ao 



0.54 



1.0 



F 

a;^ 

F 

a;^ 



(8) 



The differential conductance measured in the experi- 
ments is not exactly proportional to the density of states 
n{E) in the superconducting wire. Two effects must be 
taken into account in order to calculate dl /dV{V) from 
n{E) : Coulomb blockade and the temperature of the 
probe electrode. Coulomb blockade results from the fi- 
nite impedance of the electromagnetic environment of the 
tunnel junction The characteristics of the environ- 
ment are found from the dl /dV{V) characteristic of the 
circuit in the normal state, reached at 73 > 0.1 T, which 
presents a 10% logarithmic dip at zero voltage. The en- 
vironment can be modeled by a capacitance C = 8 fF in 
parallel with a resistance 7? = 250 fJ. Coulomb block- 
ade results in a convolution of the density of states with 
a function P{E), the probability for the electromagnetic 
environment of the tunnel junction to absorb an energy 
E El: 



§iV)^l-faEniE)PieV-E). 



(9) 
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Here, P{E) = a/ Eo{E/ Eq)"-'^ for E smaller than £^0 = 
e^/naC, with a = 2R/{h/e'^). The tunnel resistance of 
the junction was Rt — 140 kfi. As a result of this cor- 
rection, the peak value of n{E) is reduced by a few % 
in dl /dV{V). Finite temperature in the normal probe 
results in a further convolution with the derivative of a 
Fermi function. In our experimental setup, this temper- 
ature is slightly voltage-dependent, because the probe 
electrode is thermally isolated from the larger contact 
pads by superconducting connections. Heat transport 
only occurs by electron-phonon coupling and by electron 
tunneling through the junction. Since both mechanisms 
are very inefficient, even an input power Vin in the fW 
range can induce a significant temperature increase. At 
bias voltages large compared to the superconducting gap, 
heating by the tunneling current has a sizeable effect. 
In contrast, at bias voltages V slightly below A/e, only 
quasiparticules at energies larger than A — eV can tun- 
nel, resulting in evaporative cooling 0|. The effective 
electron temperature T is obtained by solving the heat 
equation: 

mT'-T^^,)-Vin+ J^E-^ n{E+eV){l-f{E)) = 0. 

(10) 

The first term describes heat transfer to the phonon bath, 
with E ~ 2 nW./m-3.K-5 for Cu [l^, =i 0.08 nw? the 
volume of the normal region of the probe electrode, and 
Tph = 25 mK the phonon temperature. The second term 
accounts for additional uncontrolled heat flow, that we 
attribute to spurious electromagnetic noise. The third 
term accounts for heat transfer throught the junction, 
with f{E) the Fermi function at temperature T. From 
the fit of the data at i? = and Is = 0, we find V-m = 
185 aW, corresponding to T = 65 mK at eV ^ Aq. 
The maximum cooling effect is reached at eV^/Ao = 0.99, 
where T = 30 mK; heating dominates for eV^/Ao > 1.02, 
with T = 210 mK at eV/Ao = 1.5. 

In Fig. 3, we present with solid lines best fits of the 
data, taking into account both Coulomb blockade and 
temperature corrections. The values of the fit parameter 
r for each curve are given in the insets. For a comparison, 
are also shown with dashed lines fits with a constant elec- 
tron temperature (T = 60 mK). The y— dependent tem- 
perature correction only matters for the sharpest curves. 
In turn, by fitting r(/s, B)/Ao with Eq. 10 and 1©, we 
find It = 240 //A and Br = 105 mT. The theoretical val- 
ues, assuming that the electrical dimensions of the wire 
are identical to the geometrical ones, are Ir = 310 /iA 
and -Br — 105 mT. Conversely, the experimental values 
of Ir oc ^q"^ and Br oc (w^o)^^ can be used to extract 
effective values ^oeS = 162 nm (instead of 125 nm) and 
WeS = 93 nm (instead of 120 nm). This corresponds 
in turn to an increased value of the diffusive coefficient: 
_D = 81 cm^s~^ and, through the resistance, to an effec- 
tive thickness ^cff = 31 nm (instead of 40 nm). Reduced 



effective dimensions for electrical transport could be at- 
tributed partly to the surface oxidation of the aluminum, 
which was exposed to air at atmospheric pressure before 
measurement, and to surface roughness. 

A by-product of the Usadel equations is a straight- 
forward calculation of the critical current. Accord- 
ing to Eq. O, /s oc Usir)dip/dx. Since ;7^(r) de- 
creases with r, Is presents a maximum as a function 
of dip/dx, which is the thermodynamic critical current. 
At i? = and fcsT ^ Ao, the maximum occurs at 
^odif/dx ~ 0.69, and corresponds, in agreement with 
[l|, to Ic « 0.75SAl^^y'N{0)a/h « 0.53/r = 125 /lA 
(using the experimental determination of Ir). The dif- 
ference with the measured Ic = 106 /iA might be due to 
the uncontrolled environment of the wire and to inhomo- 
geneities in the wire cross-section. 

In conclusion, we have measured by tunneling spec- 
troscopy on a superconducting wire the effect on the su- 
perconducting order of a supercurrent Is and of an ex- 
ternal magnetic field B. As predicted by the theory of 
superconductivity in diffusive conductors, the overall ef- 
fect solely depends on a single parameter, the depairing 
energy F. For our narrow wire, the Usadel equations 
lead to a simple expression for this depairing energy as 
a function of Is and B, which compares well with the 
experimental determination of F. 
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